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Abstract 

Recently Hofava proposed a renormalizable gravity theory in four dimensions which reduces to 
Einstein gravity with a non-vanishing cosmological constant in IR but with improved UV behav- 
iors. Here, I study an IR modification which breaks "softly" the detailed balance condition in 
Hofava model and allows the asymptotically flat limit as well. I obtain the black hole and cosmo- 
logical solutions for "arbitrary" cosmological constant that represent the analogs of the standard 
Schwarzschild-(A)dS solutions which can be asymptotically (A)dS as well as flat and I discuss their 
thermodynamical properties. I also obtain solutions for FRW metric with an arbitrary cosmolog- 
ical constant. I study its implication to the dark energy and find that it seems to be consistent 
with current observational data. 
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Recently Hofava proposed a renormalizable gravity theory in four dimensions which re- 
duces to Einstein gravity with a non-vanishing cosmoiogicai constant in IR but with im- 
proved UV behaviors 111, [J. Since then various aspects and solutions have been studied 

E SB B B B B M, H 0, H H M, M, E M, m M, B H H H Hi, H H H H2 . 

In [7J, it has been pointed out that the black hole solution in the Hofava model does not 
recover the usual Schwarzschild-AdS black hole even though the general relativity is recov- 
ered in IR at the action level. (For a generalization to topological black holes, see 



On the other hand, in [20y an IR modification which allows the flat Minkowski vacuum has 
been studied by introducing a term proportional to the Ricci scalar of the three-geometry 
jA^Ry', while considering "vanish ing " cosmoiogicai constant (~ Aw) in the Hofava gravity. 



(For related discussions, see also [lOL Il7|.) 



In this paper, I consider the black hole and cosmoiogicai solutions in the generalized model 
with the IR modification term /j, 4 R^ but with an "arbitrary" cosmoiogicai constant in the 
Hofava gravity. These solutions represent the analogs of the standard Schwarzschild-(A)dS 
solutions which have been absent in the original Hofava model. I discuss their thermody- 
namical properties also. 

To this ends, I start by considering the ADM decomposition of the metric 



ds 2 = -N 2 c 2 dt 2 + gij (dx* + N l dt) (dx j + N j dt 



(1) 



and the IR-modified Hofava action which reads 
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where 



is the extrinsic curvature, 



(3) 



(4) 



is the Cotton tensor, k, A, v, /x, Aw, and lu are constant parameters. The last term, which 
has been introduced in |2j, |l0|, |20|], represents a "soft" violation of the "detailed balance" 



condition in [2j and this modifies the IR behaviors 1 . 

Let us consider now a static, spherically symmetric solution with the metric ansatz 



ds* 



,2, 2, ,2 +r 2( de * 



-N{ryc l dt z + 



fir) 



sm 



(5) 



1 In 20], lo = 8/i 2 (3A — 1)/k 2 has been considered for the AdS case, but lo may be considered as an 
independent parameter, more generally. 
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By substituting the metric ansatz into the action (|2J), the resulting reduced Lagrangian, 
after angular integration, is given by 



if -I? 
(2A-1) U 9 -2A 



N 

8(l-3A)v7 
-2(u,-A w )(l-/-rf')-3A^r 2 



(6) 



where the prime (') denotes the derivative with respect to r. In 13], only the asymptotically 
Minkowski solution with Aw — > limit was considered. Here, I obtain the general solution 
with an "arbitrary" A^y. 

The equations of motions are 



if - l) 2 
(2A-1) U 2 -2A 



— /' + K^f 2 ~ 2(« - A*0(1 - / - r/0 - SA^r 2 = 



(£)' ((A - 1)/' - ^ + 2<c - A,),) + (A - 1)£ (r - ^) = (7) 

by varying the functions N and /, respectively. 

For the A = 1 case, which reduces to the standard Einstein-Hilbert action in the IR limit, 
I obtain 



N 2 = f = 1 + (uj - A w )r 2 - yJr[u}{u-2A w )r 3 + P] 



(8) 



where (3 is an integration constant 2 . It is easy to see that this reduces to Lii, Mei, and 
Pope (LMP)'s AdS black hole solution in [7] (I consider A^ 2 = / always, from now on), by 
identifying (3 = — a 2 /Aw, 



f = 1 - A w r 2 



a 



for uj = 0, Kehagias and Sfetsos's asymptotically flat solution in 2(|, by identifying (3 
4luM, 



(9) 



/ = 1 + cur 2 - Jr[u 2 r 3 + AuM] 



(10) 



for A w = 0. 

For r 3> [(3/uj{uj — 2A W )] 1 ^ 3 (by considering asymptotically AdS case of A w < with 
u > 0, (3 > 0, for the moment) <M} behaves as 



A 2 

/ = 1 + ^K r 2 - 
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+ 0{r~ 4 ). 



2w' 2^(0; - 2A W ) r 
This agrees with the usual Schwarzschild black hole (by adopting the units of G — c = 1 



(11) 



/=1- 



2M 
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2 If one add another IR modification term k 2 /^ 2 (8(3A — I)) 1 (3A W as in the solution becomes 

TV 2 = / = I + (uj - A w )r 2 - J r[{uj(u - 2A W ) + /3A w /3}r 3 + 0]. But this can be obtained by redefining 



the parameters Aw — > y 1 — /3/3 Aw, w — > w + (y I — /3/3 — l)Ayp in © 



for Aw = and with (5 = AuiM, independently of lu. But for Aw ^ 0, there are corrections 
in the numerical factors due to u effect: With (3 = 4uM, (iTTj) can be re- written as 



/ = 1 + M 



A 
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, 2M 1 

r 



1 + 2|Aw/o;| 



0(r" 4 ) (13) 



in which the coefficients slightly disagree with those of the standard Schwarzschild-AdS black 
hole 3 , by the factor 1 \Kw/uj\\ 

f = 1 + \*z±_™ (14) 



This solution has a curvature singularity with the power of r 3 / 2 at 



r 



^^- ^"- a 'U l), (15) 

but no curvature singularity at r = oo. Note that this singularity is milder than that of 
Einstein gravity which has R^ va ^ R^vafi ~ r 6 at r = 0. 4 

For asymptotically AdS, i.e., Aw < (with uj > 0), the solution (jSJ) has two horizons 
generally and the temperature for the outer horizon r + is given by 5 

_ 3A^r 4 + 2(o; - A w )r* - 1 

87rr + (l + (w - A w )r%) ' 1 J 

In Fig.l, the temperature T vs. the horizon radius r + is plotted and this shows that asymp- 
totically, i.e., for large r+, the temperature interpolates between the AdS cases (above two 
curves) and flat (bottom curve). There exists an extremal black hole limit of the vanishing 
temperature where the inner horizon r_ meets with the outer horizon r + at 6 



-{u -A w ) + ^{oo- A w ) 2 + 3A^ 
\l 3A£ (17) 



3 This seems to be a quite generic behavior of the broken "detailed balance" . See, for example, 

4 I consider the four-dimensional curvature invariants just for a formal reason, i.e., the comparison with those 
of Einstein gravity. But, the degree of singularity is unchanged even if the three-dimensional curvature 
invariants are considered only. 

5 Due to the lack of Lorentz invariance in UV, the very meaning of the horizons and Hawking temperature 
would be changed from the conventional ones. The light cones would differ for different wavelengths and 
so different particles with different dispersion relations would see different Hawking temperature Th and 
entropies, the Hawking spectrum would not be thermal. But from the recovered Lorentz invariance in IR 
(with A = 1), the usual meaning of the horizons and T as the Hawking temperature would be "emerged" 
for long wavelengths. The calculation and meaning of the temperature should be understood in this 
context. 

6 For Aw ^ limit, l|17p becomes 0/0. But, from (fT6|) , one can get easily = 1/\/2lo without the 
ambiguity. 
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FIG. 1: Plots of 47rT vs. black hole's outer horizon radius r + for various A\y, i.e., A\y = — 1, 0.5, 
(top to bottom) with to = 2. For large black holes, the temperature interpolates between the AdS 
cases (above two curves) and flat (bottom curve). In IR region, i.e., small black holes, there exists 
an extremal black hole limit of the vanishing temperature at in which r_ meets with the outer 
horizon r + . The region smaller than the extremal radius, which being the Cauchy horizon, shows 
an unphysical negative temperature. 

and the integration constant 

= I + 2(uj - A w y + + A 2 w rj 
r+ 

gets the minimum. The extremal radius r* + is the Cauchy horizon and so continuation to 
region of r + < r* + does not make sense to outside observer; T < for r + < r* + reflects a 
pathology of the region. (For some recent related discussions, see j30j].) 

For asymptotically dS, ie., A w > 0, the action is given by an analytic continuation 7 

H — > ifi, v 2 — > —iv 2 , uj — > — oj (19) 

of (EJ) [7]. This can be easily seen in the expansion ( |T3l) for r ^> [{3 /\uj{uj — 2Av^) |] 

Aw 



to 



r * _ _ 2M ; A + (r- 4 ) (20) 



'l + 2\A w /u 

which agrees with the usual Schwarzschild-dS cosmological solution 



/ = l-^-^, (2D 



7 This corresponds to the analytic continuation of the three-dimensional Euclidian action Wkuc = 
' Tr (r A dT + |r A r A r) + uj d 3 x^/g(R^ - 2A W ) into iW Lor with the real-valued action W Lol - 
This prescription agrees with [3l| but disagrees with [32 1. 
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FIG. 2: Plots of 4irT vs. black hole horizon radius r + for the dS black hole case with kyy = +1 
and u) = —2. There exists also an extremal limit of the vanishing temperature at of (|17p in 
which the black hole horizon r + coincides with the cosmological horizon r ++ , i.e., the Nariai limit. 
There is an infinite discontinuity of temperature at f+ = 1/y/Aw — w. 



up to some numerical factor corrections. The dS solution also has two horizons generally; 
the larger one r ++ for the cosmological horizon and the smaller one r + for the black hole 
horizon. 

The black hole temperature is given by ([Pol) also but now with Ayy > 0, uj < 0. There 
exists also an extremal limit of the vanishing temperature at r* + of (|17p in which the black 
hole horizon r + coincides with the cosmological horizon r ++ , i.e., the Nariai limit. But a 
peculiar thing is that there is an infinite discontinuity of temperature at 

r + = ~7^— (22) 
V Avk — to 

(see Fig. 2). This may be understood from the following facts. First, by writing (3 = 4ojM 
with mass parameter M, in conformity with the usual convention of (j2ip . one needs to 
consider additionally the condition 

m < ( 2A y") r 3 (23 ) 

in order that the black hole horizon exists and the curvature singularity at r = is not 
naked. But this inequality is satisfied always from the relation 

M = l + 2(u-A w y + + A^X . 
4u;r + 

( 1231) reduces to the condition [(uj — kw) r % + I] 2 > 0, where the equality for r + = 
1/^/Aw — uj = r + corresponds to the upper bound of the mass (|23l) . In other words, 
M < ( 2A ^ -a; ) r 3 f or a ^ r+ exce pt f or r + — r + , where M meets the upper bound 

Mbound " — I — r+ " 4(A w -cuy/i- (25) 



6 



FIG. 3: Plots of mass spectrum M of (|24p (lower curve) vs. the upper mass bound M = 2 ™~ r+ 
of (|23p (upper curve) for the existence of black hole horizon r + . This shows that the mass bound is 
always satisfied for all horizon radius r + and this is saturated at the point r + = l/y/A\y — to = f + . 
And also, the mass parameter gets the maximum value in the spectrum at the Nariai limit rV 
Here I plotted Aw = +1 and u> = — 2 cases but the results are generally valid for arbitrary values 
of A\y > and u < 0. 



(See Fig. 3 for the graphical explanation of this circumstance.) So, the occurrence of the 
infinite temperature discontinuity and even the negative temperature for r + < f + would 
be a reflection of being the upper bound of the mass parameter M, for a given r + , like 
the spin system with the upper bound of the energy level. (For some recent discussions in 
different contexts, see also [33].) There is no "geometrical" reason to exclude r + < r + , where 
T < 0. But, the negative temperature might be a signal of the instability of the smaller 
black hole, like the negative temperature spin systems in the ordinary surroundings with a 
positive temperature. This situation is quite different from the asymptotically AdS or flat 
case, where T < region is geometrically protected by the Cauchy horizon at r* + in which 
the inner and outer horizon coincides and T = 0. 

So far, I have studied the black hole and cosmo logical solutions for A = 1, which matches 
exactly with the Einstein-Hilbert action in IR. It would be interesting to find the more 
general solutions for arbitrary values of A. Especially in cosmology, the arbitrary A solutions 
would be also quite important for the practical purpose [20]. So I consider a homogeneous 
and isotropic cosmological solution to the action (J2J) with the standard FRW form (by 
recovering "c") 



ds 1 



-c 2 dt 2 + a 2 (t) 



dr 2 



kr 2 /Rl 



+ r 



(de 



2 i -2 

+ sin 



(26) 



where k = +1,0,-1 correspond to a closed, flat, and open universe, respectively, and 
i?o is the radius of curvature of the universe in the current epoch. Assuming the matter 
contribution to be of the form of a perfect fluid with the energy density p and pressure p, I 
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find that 



a 
a 



6(3A- r 
6(3A- 1; 
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+ 
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A 2 W 



(27) 
(28) 



I have corrected some typos in 20].) Here I have considered the analytic continuation 



for the dS case, i.e., A w > and the upper (lower) sign denotes the AdS (dS) 
case. Note that the 1/a 4 term, which is the contribution from the higher-derivative terms in 
the action (j2j), exists only for k ^ and become dominant for small a(t), implying that the 
cosmological solutions of general relativity are recovered at large scales. The first Friedman 
equation fl27j) generalizes those of [7f and [2(| to the case with an arbitrary cosmological 
constant and the soft IR modification term in 0, 10, 20]. However, it is interesting to note 
that there is no contribution from the soft IR modification to the second Friedman equation 
([28]) and this is identical to that of [7|. 

For vacuum solutions with p = p = 0, I have 



K 4 /J, 2 



16(3A - 1) 



k 



Rla 2 



-A 
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+ 



2ku 
Rfa 2 



(29) 



Here, the role of the u term is crucial. Without that term, only the constant solution of 
a 2 = — 1/K w Rq with k = — 1 exists when Aw < 0, otherwise (a/a) 2 becomes negative [7|. 
But now with the last term I have more possibilities: There may exist non-constant solutions 
even for k — +1 if —u is big enough to make (a /a) 2 > 0. Actually, for — uj > 2\A W \, k — +1, 
there exists a cyclic universe solution 



*AdS 



k(-u-\A w \) 



f?2 A 2 



1 ± 



\ 



-uj(-uj - 2\A W \) 



. ( K 2 /iA W 

sm —— r 

\2(3A - 1) 



(t-7) 



(30) 



which is oscillating between the inner and outer bouncing scale factors 



yf^2kuj ± J2k(-uj - 2\A W \) 



(31) 



2-Ro|Awi 

and the integration constant 7, depending on the initial conditions. The two bouncing scale 

factors merge as —to becomes smaller and coincide at a^ dS = yj k/\Aw\Ro when —to = 2\Aw\- 
For k = —1, the solution reduces to LMP's constant solution when u = 0, but there is no 
solution for other values when u < 0. 

On the other hand, for the dS case, i.e., A w > and u > 0, the general solution is given 

by 



Ods(*) 



2|3A-1| ± >l 
-e 2 



h 2 \ii\RqA w 
For k — — 1, a bounce occurs at 

a dS = 



(t- 7 ) , k 2 K 2 \fi\u( y uj-2A w ) T ^H^ (t „ 7) _ k(u-A w ) 



8 1 3 A — 1 |i?gAiy 



p2 A2 
n IV W 
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v /z 2EJ ± J-2k(uj - 2A 
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2R A 
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(33) 
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when u > 2A\y', at a~[g when a(t) shrinks toward a^ s , at a dS when a(t) expands toward 
aj g . uj = 2A W is the marginal case where the two bouncing scale factors coincide at a dS = 



—k/AwRo and the universe evolves monotonically from that point to de Sitter vacuum 
asymptotically or vice versa. When u < 2Aw (and also for arbitrary values of uj > when 
k = +1), the universe evolves from the big bang singularity to de Sitter vacuum or vice 
versa. For uj = 0, k = +1, this reduces to the LMP's solution with the minimum scale 
factor a min = l/y/A w R 0] 8 . 

In addition to the evolution of the universe, there are very strong constraints on the 
equation of state parameters for the constituents of our universe. So, it would be an impor- 
tant test of our Hofava gravity, whose additional contributions to the Friedman equation 
may not be distinguishable from those of dark energy, to see if one can meet the correct 
observational constraints 9 . Then, it is easy to see that the energy density and pressure of 
the dark energy part are given by (for a related discussion with matters in the context of 



the original Hofava gravity, see |34j ) 



3k 2 /j 2 ( —k 2 2kuj 2 \ 

PD - R = ± 8(3A^1) [W ~ W 2 ~ W ) ' ( } 

K 2 fl 2 ( k 2 2kuj 2 \ 

PD - E - = T 8(3X-1) [~W ~ W ~ W ) ' (35) 

respectively and the equation of state parameter is given by 

Pd.e. ( k 2 - 2kujR 2 a 2 - 3A 2 w Ry \ 
WB - E - pd.e. {3k 2 + 6kujR 2 a 2 + 3A 2 w R^J ' 1 J 

This interpolates from wd.e. — 1/3 in the UV limit to wd.e. = — 1 in the IR limit but the 
detailed evolution patten in between them depends on the parameters k,uj,A w (See Fig.4- 
7). This looks to be consistent with current observational constraints but it seems to be 
still too early to decide what the right one is [35j. But if I consider the transition point 



from deceleration phase to acceleration phase, which is given by ax = y\k\/\A w \Rl from 
(I28p by neglecting the matter contributions, the formula gives tOrxE. = —1/3, independently 
of the parameters k,u,A w . If I use ~ 1/1.03 ~ 0.9709 which corresponds to zt ~ 0.30 
in the astronomer's parametrization z = 1/a — 1, I get |A^| ~ (1.03) 2 Rj^ ~ 1.0609i?Q 2 
for the non-flat universe with \k\ = 1. And also, if I use flk ~ —0.026 [36[] in the current 
epoch (a = 1) for the deviation from the critical density, Qk = fi 2 k\Aw\Lp/2a 2 H 2 RlMp, 
Hubble parameter H = a/a, the ratio of Planck mass and length Mp/Lp = 2(3A — 1)/k 2 , 
k = —1, I get fi ~ 0.2214HqRqMp/ L P with the current value of Hubble parameter H . 10 
Finally, if I use u>d.e. ~ —1.08 in the current epoch, I get uj ~ 1.0067-Rq 2 which predicts the 



10 



8 It is interesting to note that these various scenarios for dS case have been considered earlier by Calcagni 
in the original Hofava model [3| but it is important to note that these can be "realized" only in our 
modified model with the uj terms. In Calcagni's notation, the scenarios are categorized by the values of 
3c 2 K 2 — 4£? 2 |A| but one can have only 3c 2 K 2 — 4£> 2 |A| = if one use the actual values of the parameters 
for the Hofava model, i.e., (34) and (35). This corresponds to the LMP's solution. 

9 While this paper was being finalized, I met a paper by Mukohyama [28 1 which propose the dark "matter" 
as integration constant in Hofava gravity. 

I follow the physical convention of Ryden 37| which disagrees with [l|, |2j . 
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FIG. 4: Plot of equation of state parameter u>d.e. vs - scale factor a(t) for u> 2 > A 2 ^, kuj < 0. 



There are two infinite discontinuities of wd.e. at = y — ku> ± |A;| yuj 2 — A^/|Avk|-Ro where pn.E. 
vanishes. Here, I considered \u\Rq = 2,\A\y\R 2 l = 1 case {ojRq = —2,k = +1 or ujRq = +2,k = 
-1). 

evolution of wd.e. as one of the curves in Fig. 6 since uj < |Aw|: IfIusei?o ~ 6.2017 c/H from 
fi fc = kc 2 /H 2 R 2 0.026 and H ~ 70km s^^pc" 1 , I get A w ~ 1.5018 x lO^Mpc" 2 , u ~ 



1.4251 x 10- 9 Mpc" 2 , fi ~ 5.6636 x 10 35 kg s -1 . 

Note added: After finishing this paper, a related paper (3§| appeared whose classification 
of all the possible cosmology solutions in the Hofava gravity without the detailed balance is 



overlapping with mine. (See also [39].) 
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